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We investigate experimentally the effects of light assisted collisions on the coherence between mo- 
mentum states in Bose-Einstein condensates. The onset of superradiant Rayleigh scattering serves 
as a sensitive monitor for matter wave coherence. A subtle interplay of binary and collective effects 
leads to a profound asymmetry between the two sides of the atomic resonance and provides far bigger 
coherence loss rates for a condensate bathed in blue detuned light than previously estimated. We 
present a simplified quantitative model containing the essential physics to explain our experimental 
data and point at a new experimental route to study strongly coupled light matter systems. 
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Elastic Rayleigh scattering of photons from atoms in a 
Bose-Einstein condensate (BEC) creates long-lived rip- 
ples in the density distribution of the atomic cloud. 
Bosonic stimulation leads to a positive feedback mech- 
anism enhancing the formation of a matter-wave grat- 
ing which scatters photons coherently predominantly 
along the directions of high optical depth. This directed 
Rayleigh scattering is well known as Rayleigh superradi- 
ance and has been studied in BECs extensively in differ- 
ent geometries [IH3] ■ The coupled dynamics of superradi- 
ant (SR) scattering with simultaneous build-up of recoil- 
ing matter-waves and light fields have been successfully 
described using Maxwell-Schrodinger equations as well as 
rate equations derived from those [5]. 

Recently, those established descriptions of SR scatter- 
ing have been challenged by a new model predicting a 
peculiar asymmetry in the dependence of the dynamics 
on detuning of the drive light [5] . Clear experimental evi- 
dence of such an asymmetry, which is pronounced at high 
atomic density and vanishes at low density, is reported in 
[7] . This has sparked an ongoing debate [8] and it appears 
that the mechanism leading to the different dynamics for 
below resonance (red) and above resonance (blue) tuning 
of the incident light is not yet fully understood. 

To address this open question we present here an ex- 
perimental study of the threshold behavior of SR light 
scattering for a wide range of detunings up to 35 GHz. 
The observed detuning dependence rules out previous ex- 
planation attempts based on the action of dipole forces 
[9] . We offer a physically motivated explanation for 
the asymmetry in the threshold behavior of SR scatter- 
ing based on detuning dependent loss of matter-wave co- 
herence. Resonant excitation of close pairs of atoms to 
excited state molecular potentials and subsequent spon- 
taneous decay provides a source of frequency shifted pho- 
tons, which for blue detuned drive light can be trapped 
inside the BEC for a long time. The different accessible 
molecular branches with continuous and discrete spec- 
tra together with the fact that the inelastically scattered 



photons are inevitably red shifted with respect to the 
excitation light lead to an intrinsic red-blue asymmetry 
in the loss rates. We compare the results of a simplified 
rate equation model including such loss to our experimen- 
tal data and find satisfactory agreement. The underly- 
ing mechanism offers exciting new possibilities to study 
optical excitations trapped deep inside dense ultracold 
atomic gases. 

In our experiments, SR Rayleigh scattering is induced 
in a trapped BEC by illuminating it with an off-resonant 
light pulse along the long axis of the condensate, which 
leads to recoiling atoms gaining two photon momenta. 
We use a 100 /is rectangle pulse and vary the pulse in- 
tensity to explore the onset of SR. BEC's of 87 Kb atoms 
in the \F — 1, mp — —1) hyperfine state are prepared by 
evaporative cooling in a Ioffe-Pritchard magnetic trap 
[3]. We obtain prolate condensates containing 2.2 x 10 5 
atoms with in-trap Thomas-Fermi radii of rj_ rj 6 /im 
and ry ~ 60 /im in the radial and axial directions, with 
no discernible thermal fraction. The pump light is de- 
tuned from the \F = l,m F = -1) \F' = 2,m F , = -2) 
transition on the Di line of 87 Kb at 795 nm and is 
circularly polarized. We measure the drive detuning 
(A = -2.57 35 GHz) using a wavemeter (~ 300 MHz 
resolution) and a Fabry-Perot resonator (1.5 GHz FSR) 
referenced to a laser stabilized to saturated absorption 
features of the 85 Rb Z?i line. The pump beam is focused 
to a waist radius of 20 fim (at e~ 2 ) on the atoms. The 
pump pulse duration is chosen long enough to suppress 
backward (Kapitza-Dirac) scattering of atoms and short 
enough to neglect decoherence due to decaying overlap 
of matter wave packets in the later theoretical modeling. 
The range of explored pump detunings is chosen such 
as to have negligible pump light depletion at low detun- 
ings and is limited by our available laser power at high 
detunings. 

The populations of atomic momentum modes (Ohk & 
2hk) are extracted from absorption images taken after 
time-of- flight. At each detuning we measure the popula- 
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FIG. 1. (color online) The transfer efficiency as a func- 
tion of normalized single atom Rayleigh scattering rate at 
A = 5.24 GHz; Open triangles: population of the Ohk mo- 
mentum mode; Filled circles: population of the 2hk momen- 
tum mode; Dashed line: linear fit; Inset: time of flight image 
of the atomic momentum distribution. 
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FIG. 2. (Color online) Normalized threshold scattering rate 
vs. detuning. Filled symbols show the measured threshold 
increase; Error bars designate 95% confidence level; Gray 
shaded region depicts the expected threshold due to close 
range dipole-dipole interaction allowing for a factor of two 
variation of the light induced loss rate L ge 



tion transfer for different single atom Rayleigh scattering 
rates R, which in simple models determine the timescale 
for the dynamics. We fit a straight line to the results 
with low scattering rate to extract a phenomenological 
threshold pump rate i?t^,(A) where superradiant gain ex- 
ceeds linear losses enough to start significant population 
transfer during the interaction time. Figure [T] shows an 
example of transfer efficiency measurements at a blue de- 
tuning A = 5.24 GHz. The scattering rate is normalized 
to the measured threshold rate Rq at a red detuning of 
A = —2.57 GHz [10]. A detuning asymmetry in the 
threshold and a saturation of the transfer efficiency at 
high scattering rates is visible. 

In Fig.[2]we present measured threshold rates as a func- 
tion of pump laser detuning. A threshold increase up to a 
factor of three is evident for low values of the blue detun- 
ing. Also shown in the figure is the expected threshold 
increase due to light assisted collisions calculated from 
the model presented in the following. 

The starting point for the description of the coherent 
part of the SR process are coupled Maxwell-Schrodinger 
equations From these, we derive a rate equation for 
the number of atoms N 2 appearing in the recoil mode 
N 2 = G ■ (N 2 + 1), to describe the early stages of the 
dynamics where depletion of condensate atoms Nq and 
pump light is not important [9 . The rate constant for 
growth G = Rbo depends for a fixed sample geometry lin- 
early on the single atom Rayleigh scattering rate R and 
on the effective resonant optical depth & of the sam- 
ple along the propagation direction of the superradiant 
light mode [5]. Since the Rayleigh scattering rate varies 



symmetrically with laser detuning no asymmetry is pre- 
dicted by this model [12] . At this level of description 
Rayleigh superradiance does not have a threshold pump 
rate. To have a more realistic model of the onset, damp- 
ing mechanisms for the coherence gratings need to be 
accounted for. We include loss rates Lr = R to describe 
removal of N 2 atoms by spontaneous Rayleigh scatter- 
ing, L gg (n$) to describe damping by incoherent ground 
state collisions between (N ,N 2 ) pairs [T3] , and a loss 
rate L ge to account for light assisted collisions followed 
by radiation trapping. We also add a nonlinear loss rate 
L n i(No, A), without discussing here the physical origin 
of possible nonlinear loss terms further [9J , to account 
for processes that depend nonlinearly on the population 
of the superradiant modes as those suggested in [7] . The 
resulting rate equation 

N 2 = G{N 2 + 1) - (L R + L gg + L ge )N 2 - L nl N% (1) 

shows initial exponential growth when gain exceeds lin- 
ear losses, i.e. G > Lr + L gg + L ge . Equality of gain 
and linear losses defines the threshold for SR scattering 
in this model. The threshold does not depend on L„;. 
Similar to a depletion term, L n i clamps the growth rate 
later during the evolution when the population in the re- 
coil mode becomes significant. We denote the threshold 
gain (Rayleigh rate) in the absence of the L ge term as 
Go {Rq) and parametrize the light induced loss rate as 
L ge = Rxno(l + Fn). Here, i?x«o is the light assisted 
collision rate, with n the condensate density and x a 
molecular parameter, F is the fraction of resonant pho- 
tons produced in a collision, and n is the average number 
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of subsequent scattering events for a resonant photon in- 
side the cloud. 

This parametrization of L ge puts emphasis on the role 
of N 2 atoms for the contrast and spatial coherence of 
the matter wave grating responsible for the amplified di- 
rected SR scattering. The matter wave grating amplitude 
can decrease or dephase either by direct participation and 
subsequent loss of an N 2 atom in a binary light assisted 
collision or by interaction of an N 2 atom with a resonant 
photon produced in any light assisted collision within the 
cloud. As we will show later, this second mechanism is far 
more important than direct loss of N 2 atoms [H] . Since 
in a fully quantum mechanical picture of SR scattering 
recoiling atoms and backscattered photons are created 
as correlated pairs contributing on equal footing to the 
gain, it is interesting to ask how important the loss of 
photon coherence due to inelastic collisions is for the net 
reduction of SR gain. The vast majority of binary colli- 
sions, which are the source for nearly isotropic incoherent 
resonant radiation, happens between N atom pairs as- 
sisted by pump light photons. As long as depletion of 
the pump light is negligible, the direct influence of the 
frequency shifted radiation on the coherence of the light 
grating is marginal. It is the strong response of atoms 
to even minute amounts of resonant light which spoils 
the coherence of the matter wave and this way also the 
mutual coherence between light and matter waves. 

Using Eq. ([!]) and the above expression for the loss 
term L gei the expected change of the threshold pump 
rate in the presence of light assisted collisions can be 
expressed as 

Rth(&) = L_ xn (l + Fn) \ - 1 
R V b o ~ 1 J 

which is the quantity determined in our experiment. Due 
to the strong dependence of n on optical thickness, dis- 
cussed later, the threshold increases markedly for dense 
and optically thick clouds, while it is unaltered in the 
limit of low density and optically dilute clouds, in accor- 
dance with the experimental observations in [7J. 

To allow for a comparison between model predictions 
and our experimental data, we turn now to a more de- 
tailed discussion of the three-step process leading to the 
coherence loss rate L ge and give quantitative estimates 
of the microscopic parameters X: F an d n. Modeling 
off-resonant light scattering as purely elastic ceases to be 
a good approximation at high atom densities when the 
probability to excite close pairs of atoms becomes sig- 
nificant. To include this into the description a molecular 
point of view is necessary. Here, photons can be scattered 
with a significant frequency shift with a concomitant 
change in kinetic energy of the outgoing pair of atoms, 
high enough for both atoms to leave the trap. Since the 
outgoing photon has a frequency close to the atomic res- 
onance it will scatter repeatedly inside the cloud before 
leaving. The dependence of L ge on detuning is inherently 
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FIG. 3. (Color online) Ratio x n o of pair collision rate to 
Rayleigh scattering rate for an atomic density of 10 14 cm' 3 
as a function of detuning. Dashed and dash-dotted lines: Sep- 
arate contributions from the Of and l u molecular potentials; 
Solid line: total ratio. Inset: Typical ground and excited 
state potentials for cold collisions in light fields in the de- 
tuning regime where the excited state structure is dominated 
by resonant dipole-dipole interaction. Red detuned light ex- 
cites pairs to bound states (free-bound transition), while blue 
detuned light excites to a continuum (free-free transition). 

asymmetric since for red detuning only discrete bound 
molecular states can be resonantly excited (photoassoci- 
ation resonances) while for blue detuning a continuum of 
states on repulsive molecular potentials is accessible (ra- 
diative heating) as sketched in the inset of Fig. [3] The ini- 
tial light assisted collision step in the three-step process 
has been used in the past to assess light induced atom 
loss rates from condensates [15] . To estimate quantita- 
tively the event rate coefficient for binary light assisted 
collisions R\ appearing in L ge we switch to a microscopic 
description of the collision employing the methods out- 
lined in |16j . Using a reflection approximation the event 
rate coefficient R\ is written as: 

RX = x 4tt 2 K 2 x -L \* g (R c ,E)\ 2 (3) 
Moo F>c 

Here, \i is the reduced mass for a 87 Rb atom pair, 
hkoo is the relative momentum in the entrance chan- 
nel with corresponding kinetic energy E, and R c is 
the Condon radius where the molecule tunes into res- 
onance. The Franck-Condon factor \V g (R c ,E)\ z /D c , 
with ty g (R c ,E) the ground state scattering wave func- 
tion and D c the difference in potential slopes, regulates 
the detuning dependence. The radiative coupling poten- 
tial V c = V eg (R n ) = bc{R n )fi^A, where £Ia is the atomic 
Rabi frequency and bc(R n ) a molecular parameter re- 
flecting the change of the electronic wavcfunction with 
internuclear distance R n , varies only little in the range of 
atomic distances relevant here. The ground state scatter- 
ing wave function is calculated numerically by integrating 
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the Milne equation. The accessible repulsive excited state 
molecular potentials are parameterized with dispersion 
coefficients 6*3(0+) = 11.9 a.u. and Cs(l u ) — 5.89 a.u. 
following [T7]- The main panel of Fig. [3] shows the ra- 
tio x n o between light assisted collision rate and isolated 
atom Rayleigh scattering rate at a typical BEC density of 
iiq = 10 14 cm~ 3 for excitation to the different molecular 
potentials. 

To model the second step, the frequency shift of pump 
photons assisting a collision, we estimate the frequency 
spectrum of light exiting the collision complex and deter- 
mine the fraction F emitted within one natural linewidth 
of the bare atomic resonance. Using a semiclassical wave 
packet approach we calculate the trajectory of an ex- 
cited atom pair separating from the Condon radius and 
use energy conservation to determine the frequency shift 
of light emitted along the way [9] . Due to the high ac- 
celeration on the repulsive molecular potential almost all 
of the atom pairs reach the asymptotic kinetic energy be- 
fore decaying, making the frequency redistribution func- 
tion sharply peaked around the atomic resonance. Tak- 
ing into account hyperfme branching in the decay around 
43% of the emitted light is fully resonant with the hosting 
cloud for all A > used in the experiment. 

As the third ingredient to assess the damping of the 
matter wave coherence we need to find the average num- 
ber of scattering events n for resonant photons before 
leaving the cloud. Neglecting further frequency redis- 
tribution we use a simplified Holstein model, essentially 
a diffusive transport equation for light intensity in a 
medium of high optical depth, to describe radiation trap- 
ping [18] . We use the decay time of the slowest Holstein 
mode Tq 1 = "fb 2 T nat ~ fiT nat to calculate n [19] . Here, 7 
is a geometry parameter, r nat is the radiative lifetime of 
the excited atomic state, and b is the optical depth. For 
a simple estimate we assume a Gaussian spherical geom- 
etry (7 ~ 0.06) with an optical depth equivalent to the 
geometric mean along the different condensate axes. For 
our parameters we find n > 1000. 

The gray area in Fig. [2] depicts the calculated thresh- 
old increase via Eq. (J2| for the experiment. Given the 
several rather crude assumptions in the calculation of 
the microscopic model parameters, together with smaller 
systematic uncertainties in the experimental parameters 
atom number and density, we allow for a factor of 2 vari- 
ation of the calculated loss rate L ge in Fig. [2] We note 
that the exact spectral signature is sensitive to molecu- 
lar hyperfme structure, which is not taken into account 
in our model potentials. While the qualitative and near 
quantitative agreement between data and prediction is 
satisfying to see, the rate equation and radiation trap- 
ping model applied does not do full justice to the under- 
lying complicated many-body physics. Simple inspection 
of the threshold condition reveals that at our highest ob- 
served threshold increase more than 60% of the atoms 
should have interacted with a trapped resonant photon. 



The little observed recoil heating in the experiment is 
clearly incompatible with a picture of individual atoms 
receiving random recoil kicks from an isotropic radiation 
field. In fact, at resonance and high density photonic and 
atomic degrees of freedom mix strongly, forming polari- 
ton type excitations with an effective mass very different 
from the bare atomic mass [20] . In using the radiation 
diffusion model we implicitly assume that the dephasing 
rate is still governed by the bare atomic decay rate F. 

Whatever the precise nature of the trapped excitations 
is, their incoherent production by collisions and slow dif- 
fusion implies the presence of electronic excitation inside 
the cloud many natural lifetimes after the pump light has 
left the cloud. Photoionization out of the excited state 
can provide a critical test of the model but also a tool to 
study the temporal dynamics of the polaritons in detail. 
The spatial structure of phase damage can be mapped 
out by matter wave shearing interferometry |21j . In fu- 
ture work, it will be interesting to study the different 
response of dense clouds to resonant light applied from 
the outside or created directly inside the medium [2"2"] . 

In conclusion, we have studied experimentally the 
threshold asymmetry of Rayleigh superradiance. We 
have developed a simple, yet quantitative, model to ex- 
plain our data and discussed the underlying physics, 
which can hopefully serve as a useful guide for more rig- 
orous theoretical studies. The mechanism for deposition 
of resonant photons deep inside a dense cloud might pro- 
vide a promising route to observe Anderson localization 
of light in cold atoms [23] . 

We acknowledge support by the EU through FP6 
RTN EMALI and through the projects HIDEAS and Q- 
ESSENCE. We thank A. S. S0rensen for inspiring discus- 
sions. 
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SUPPLEMENTARY MATERIAL 
Derivation of the rate equation 

In the main text we use a simplified rate equation for the number of recoiling atoms to describe the onset of 
superradiance (SR) in the absence of incoherent losses. Here, we sketch the steps to arrive at the rate equation, which 
averages over propagation effects, starting from Maxwell-Schrodinger equations, which include all propagation effects. 
The basic coupled evolution equations for light and matter fields have been presented in the literature already several 
times 131 [11]. On the way to the rate equation we discuss, in particular, the various approximations that enter the 
derivation, also in view of recent attempts to explain the detuning asymmetry in superradiance based on approximate 
analytic solutions of the Maxwell-Schrodinger equations and modifications thereof [51 [7] . 

The starting point are the mean-field Gross-Pitaevski equation for the matter field ip in the electronic ground level 
including an effective coupling term to the light field and a classical wave equation for the propagating electric field 
E with a polarization term to describe the coherent radiation by the driven atoms: 

d h 2 , (d+ • E") (d" • E+) , l2 



* |>. (5) 



dV 2M r hA 

~^dt 2 J" ~ c 2 e dt 2 

Here, superscripts ± denote positive and negative frequency components, M is the atomic mass for s7 Rb, h is the 
reduced Planck constant, d is the atomic dipole matrix element on the driven electronic transition, c is the speed of 
light, and eo is the free space permittivity. The nonlinear term proportional to go describes the mean field energy 
from ground state van der Waals interaction. The macroscopic polarization is given by = — \ip\ 2 d (d • E*) / (HA). 
Several approximations have already been applied to arrive at this form of the equations. The trapping potential for 
the atoms has been dropped, since it has negligible influence on the dynamics on the time scale of the interaction 
with the pump pulse. Implicitly the effect of the trapping potential is, of course, contained in the initial density 
distribution of the cloud. Only the Rayleigh scattering channel back into the initial Zeeman sublevel is considered for 
the superradiant dynamics, since the inhomogeneous magnetic field destroys rapidly the coherence between different 
sub-levels. More importantly, the dipole response of the atoms is calculated in second order perturbation theory for 
an isolated atom in a rotating wave and low saturation approximation to eliminate excited states adiabatically. The 
radiative damping term iT/2, which formally needs to be added to the detuning A, has negligible influence for the 
range of detunings considered later on, and is thus left out. Doing the adiabatic elimination at the single atom 
level, the light mediated interaction between atoms in their near-field and its impact on the scattering properties are 
neglected. We choose to take these effects into account later by switching to a molecular picture for close pairs of 
atoms. Effective macroscopic descriptions including the near-field interaction, derived several times in the literature, 
have many subtleties |24[|25) . In particular, frequency redistribution processes and radiation trapping, invoked as the 
source of decoherence in the main text, are buried deep under the formalism and are only hard to recognize in the 
macroscopic treatment, which concentrates on the stationary linear response of the scattering medium. 

At this stage, as long as the Maxwell-Schrodinger equations are solved simultaneously, all coherent light mediated 
interaction between driven dipoles in their far-field and all dipole forces are still fully accounted for. Keeping the 
full 3-D description throughout is cumbersome, so the next step is a reduction to 1-D. Two potentially important 
effects, depending on the specific geometry of an experiment, that are lost in a 1-D description are the transverse 
components of the dipole force exerted by the light field distribution on the atoms and the corresponding back action 
on the light (diffraction and lensing). In our experiments, the transverse dipole force accelerates the atoms radially 
outward, increasingly for higher blue detunings. At the threshold for SR scattering the peak density is calculated to 
drop, during the short interaction time, by about 3% for the lowest blue detuning, while for the highest blue detunings 
a density change of about 20% is expected. This systematic effect is not taken into account in the data analysis and 
model presented in the main text, but can possibly explain why the observed threshold rates at high blue detunings 
dip below the red detuning reference value, via a reduced rate of incoherent ground state collisions. Lensing of pump 
light due to the spatial variation of the refractive index, modifies the intensity distribution and hence the effective 
light-atom coupling. This effect can mimic the observed detuning dependence of the threshold pump rate. For the 
lowest blue detuning in our experiment, where the effect is biggest, we estimate, from numerical 3-D simulations of 
light intensity inside the sample, a relative change of the Rayleigh scattering rate of 20% with respect to the red 
detuning reference value, thus significantly smaller than needed to explain the experimental observations by lensing 
alone. We have not included this systematic effect in the data presented in the main text. 
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This said, we simplify the above equations by going to an effective 1-D geometry assuming a constant transverse 
cross section A of the interaction region. This keeps longitudinal components of interaction and dipole forces only. To 
simplify the equations further we split the electric held into forward and backward propagating modes and introduce 
slowly varying envelope functions for the modes denoted by subscripts ± in the following. Likewise, the matter wave 
function is split into recoil modes with slowly varying envelopes as tp(z,t) — ^2 m =2n ^Pm(z, t)e~ l/yU}mt ~ mkz ^ where 
m = 2n and n is an integer number (SR order number), u m = m 2 uj r , co r = ||y is the recoil frequency, and k is the 
wave number. Since the condensates used in the experiments are much longer than an optical wavelength (L ~ 100A) 
the mode functions are orthogonal to a very good approximation. As the next step we transform the equations to 

dimensionless form, by defining electric held, time and length units via E± = e± y |^ • ^f, r = 1uj r t and £ = kz. 
Using the slowly varying envelope approximation the transformed equations read: 

de + 



^ = ,A]T { £+ V4+ 2 ^e 2 ^+^ + e_ IVU 2 } (7) 



dip m i d 2 ip m dip 

m • k ( \ 1 2 . i 1 2 \ / - F \ 1 \ ^ 



m- 



dr 2 dt? d£ 



zA (| £+ | 2 + M 2 ) ^ - i^Y.Y.^-i^.e-^- " < 

r n I 



-iKe* + e-e- 2 ^ m+1 ^^ m+2 - iAe*_e+e 2 ^ m ~^ i\> m _ 2 . (8) 



The coupling constant is expressed as A = Too/ (4A^4) and comf denotes the mean-held ground state interaction. 
The ratio of atomic absorption cross-section ao to sample cross section A and the detuning A in units of the linewidth 
r determine the strength of the effective atom-light interaction. The first two terms on the r.h.s. of Eq. Q describe 
wave packet spreading and recoil induced drift of the matter wave envelopes, which ultimately leads to coherence loss 
by spatial separation. In our experiments the influence of these terms is small due to the short interaction time and 
long sample length, hence we drop them in the following. 

To connect the equations more directly to observable quantities we switch to density matrix elements instead of 
mode amplitudes. Since we are interested mainly in the onset of superradiance we restrict the number of recoil 
modes for the matter wave to the first two. This makes the system formally equivalent to a coherent two-level 
amplifier/absorber, with a weak nonlinear contribution due to the mean field interaction. 

° h '"■ - iA (e*_e + e 2lT M4 - <£_e- 2 ' v ^ 2 ) (9) 



Or 



d ^ ~ ik(e*_e + e 2lT r 2 ^-£* + e-e- 2lT Ml) (10) 



dr 

9 ^2 = iAe*_e + e 2 " (|V- 2 | 2 - |^ | 2 ) + ^ (|^| 2 - |^o| 2 ) ^ (U) 



Or 



The last term in Eq. (Ill stems from the extra energy cost to create a density modulation in the interacting cloud 



and describes e.g. the mean field shift of a Bragg resonance. The light flux propagation equations read: 



file. 



2 



d\e-\ 



iA (r 2 ^ e 2tT £-£+ ~ i>l^>2e- 2ir e* + E_) (12) 
iA (r 2 ^oe 2 "e*_e + - ^ ^ 2 e~ 2iT e* + e-) (13) 
9 e*_e+ = -iAr ^e~ 2lT (\e + \ 2 + |e_| 2 ) - 2zA (Vo| 2 + |V^| 2 ) e*_e+. (14) 



The last term in Eq. ( 14 ) reflects the modification of light wavelength due to the refractive index of the cloud. Writing 
the complex coherences in polar form as ele+ = pie 1 ^ 1 and V'oV'2 — Pa,e lc ^ a brings the equations into a form suitable 
for further discussion. 

d d 

|V'2| 2 = -^|^o| 2 = 2Ap i p a sin(2r + ^-0 a ) (15) 



8t dr 
-\£ | 2 = - 



'e_| 2 = — \e+\ 2 = -2Ap a p i sin(2r + i -0 a ) (16) 
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= -A (| £+ | 2 + |e_| 2 ) Pa sin(2r + fa - tf> a ) (17) 
^p a = A (| Vol 2 - 1^2 1 2 ) Pi sin (2r + ^ - O ) (18) 



One can recognize the first equality in Eq. (151 as the local conservation of atom number, which is a consequence of 



neglecting wavepacket drift and spread. Similarly, the first equality in Eq. ( 16 ) expresses a continuity equation for 
the photon density, a necessary consequence of the adiabatic elimination of excited atomic states. 

We note that the terms describing the mean-field interaction and the refractive index drop from the magnitude and 
coherence equations. These terms will weakly influence the evolution of the phase of matter and light gratings (<j) a & 
fa). With our choice of the pump laser frequency as the carrier frequency for both forward and backward propagating 
light modes, the time dependence of fa acquires the recoil shift of the backscattered light and compensates the explicit 
time dependence in Eqs. ( 15fl8 1. 



To model the onset of superradiant scattering we make use of the specific initial and boundary conditions in the 

experiments, i.e. I^ol 2 3> IV^I 2 and \e+\ 2 S> | er | 2 . At first sight the equations seem to imply an odd symmetry in 

the dependence on the sign of pump laser detuning (A oc A" 1 ), and hence to explain the observed asymmetry in the 
SR threshold. A closer inspection reveals, that this is not the case. Due to the sinusoidal dependence on the relative 
phase of light and matter wave gratings, Eqs. ( 15f 18 1 support runaway solutions, growing nearly exponentially in time 



and in space, for both signs of the parameter A (A). The SR is triggered by spontaneous Rayleigh scattering that 
creates random gratings and provides a seed for the growth [26] , 

If boundary conditions are such that both |e_| 2 and |e+| 2 are strong light fields incident on the sample, the coupled 
equations describe just Bragg diffraction of matter waves in a (walking) standing wave including the backaction of 
atoms onto the light field. The different relative phase for gain between light and matter interference patterns for 
blue and red detuning can be understood easily in an optical lattice picture, when considering the spatial structure of 
Bloch waves at the band edge. Similarly, for initial conditions such that both atomic recoil modes arc macroscopically 
populated, reflection of light from a density grating is described. 

Returning to SR scattering, in a minimalistic approach the equations can be reduced to a zero dimensional system 
leading to the rate equation for the number of atoms appearing in the recoil mode. To do this we assume initial 
homogeneous matter wave coherence p a (r = 0) over the sample corresponding to one delocalized atom in the recoil 
mode and perfect phase matching conditions sin(2r + fa — <p a ) ~ 1 which is valid in the early stages of the dynamics. 



Now, Eq. (17) is solved subject to the boundary condition that pi(£, m ax) = 0. The result is inserted into Eq. (15 1, 
leading to: 

>i 2 = 2A2|e + |2 4-£)- (19) 

Integrating this over the length of the sample renders a rate equation for the number of atoms in the recoil mode as: 

^-N 2 =A 2 N Q \e + \ 2 (N 2 + l). (20) 

OT 

Restoring physical units the gain constant appearing on the r.h.s. can be written now as G = Rba, with the single 
atom Rayleigh rate R and the on-resonance optical depth b as stated in the main text. 

We conclude this derivation with a brief discussion of basic scaling properties of the dipole force and point in this 
context at some fundamental flaws in the model put forward by Deng et al. [7] to explain the observed detuning 
asymmetry in SR scattering. In a low saturation approximation the dipole force can be written as the gradient of 
the polarization energy of the atoms. For the corresponding potential we have, thus, Udip oc I/A. Here, / is the total 
intensity at the location of an atom. It is instructive to consider the detuning dependence of the dipole potential 
Udip for the case of constant Rayleigh scattering rate R, which is of relevance for the pump light in SR scattering 
experiments. Since R oc I/A 2 , we must have Udip oc A, i.e. the dipole force and potential increase in this case linearly 
with the detuning. Finally, for the onset of SR scattering the relevant dipole potential stems from the interference 
pattern between pump light at constant Rayleigh rate and scattered light. Here, the scattered light electric field E s 
is independent of detuning, while for the pump light field strength we have E p oc A, which means that the strength 
of the dipole force is independent of detuning for this case. This reflects nicely the fact that in both coupled wave 
and rate equation models for SR the Rayleigh rate is the only relevant pump light parameter for the dynamics and 
that the very mechanism of SR Rayleigh scattering can be understood in terms of dipole forces. From these scaling 
arguments we conclude that the experimentally observed detuning dependence of the threshold asymmetry cannot be 
reconciled with models based solely on the action of dynamically evolving dipole forces. 
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FIG. 4. (Color online) Left figure: Fraction of photons emitted in a frequency interval of width F above the atomic resonance 
as a function of detuning. Right figure: emission spectrum (4^) as a function of u — (uii — uj) /A, evaluated at A = 3GHz. 

In the model of Deng et al. a dipole potential is invoked as being produced by a seed photon pulse, that is 
amplified upon propagation through the pump light dressed medium [37]. This ignores both the intensity of the 
pump light, as well as the interference term between pump light and scattered light. Even in the most favorable 
case of plane wave pumping the dipole forces due to the interference term are orders of magnitude bigger than the 
contribution of the scattered light intensity alone. Due to this neglect Deng et al. also arrive at an incorrect scaling 
with detuning. Deng et al. suggest furthermore in their model, that the structure factor for light scattering is directly 
modified by their invoked incomplete dipole potential. Since a structure factor depends sensitively on interaction and 
correlation properties of the constituent particles, this suggestion appears entirely unfounded without a detailed look 
at microscopic properties. 



Frequency redistribution function 

To get the frequency spectrum of photons produced in binary collisions, we use the Ehrenfest theorem to calculate 
by classical mechanics the trajectory of the excited state wave packet in the repulsive molecular potential. Knowing 
the kinetic energy as a function of time along the trajectory, allows to transform the probability distribution for 
decay as a function of time to the spectral distribution. The wave packet approach is justified by the same stationary 
phase argument that is used to calculate the Franck-Condon factor for the upward transition. The calculation is 
completely analogous to the survival probability estimate used in the Gallagher-Pritchard model for binary collisions 
in red detuned light fields [28] . 

We start by considering the total energy E available in a generic two-body (half-) collision in a repulsive r -3 
potential, 

E = T + V (21) 
HA=lr* + % (22) 

with T and /j, the kinetic energy and reduced mass, respectively. We do not consider a centrifugal potential term. 
The ground state scattering wave function has s-wave symmetry, while the electronic angular momentum coupling in 
the non-centrosymmetric dipole potential is accounted for by the designation of the molecular state, parametrized by 
C3. The energy of accessible (1=0, 1=1) rotational state continua differs for the smallest Condon radii by less than 
the atomic natural line width hT which is dwarfed by the total collision energy A ^> T. The detuning is defined as 
A = ljl — wo, where wl is the laser frequency and ujq is the atomic line resonance frequency. 

Simple algebraic manipulation allow us to derive the differential equation describing the temporal change of kinetic 
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energy along the trajectory as 



T 



18 



1/2 



T^ 2 (E-T) 4 / 3 , 



(23) 



which can be integrated by separation of variables. Introducing scaled variables u = T/E and r = T m t, where 1/I\ 
is the radiative lifetime of the excited molecular state, we write the solution as: 



t(u) 



dx 



aT J x 1 /2(l - X )4/3 



(24) 



The coefficient a = 6(A/r) 5 / 6 (2a; r /r) 1 / 2 (fir/(C3/c 3 )) 1 / 3 contains all physical parameters of the specific system, while 
the integral can be expressed in terms of hypergeometric functions. When evaluated numerically, care must be taken 
to treat the singularities of the integrand correctly. 

The probability density for decay to the electronic ground state of the colliding atom pair is given by 



dP 
~dt 



T m exp(-T m t). 



(25) 



Moving to the scaled kinetic energy as the independent variable and transforming the differential accordingly we arrive 
at 



dP 



1 



du oT 1^/2(1 _ u )4/3 eX P [ aT J Q x l/2(l _ 2)4/3 



r„ 



dx 



(26) 



Recognizing that the variable u £ [0..1] maps the energy of the outgoing photon on the interval [ojq + A..wo] the 
fraction F of photons emitted in a frequency interval of width T above the atomic resonance can be written as 



F = exp 



r™ f 1 -^ dx 

aT J X 1 /2(l_ a; )4/3 



(27) 



The spectrum of fluorescence for excitation at A = 500r as well as the fraction Fasa function of detuning are shown 
in Fig. [4] The additional (small) broadening due to the finite emission time T" 1 is not taken into account in this 
simple calculation. We finish this calculation with the remark, that the recoil shift for the emitted radiation, which 
is, of course, negligibly small compared to the red shift compensating the change of relative kinetic energy, must be 
evaluated using the total mass of the composite radiating system, i.e. twice the atomic mass. Providing the answer 
to the equivalent questions about deposited energy and recoil for the case of resonant radiation incident on a whole 
group of close atoms, is an interesting but highly nontrivial task, in our view. 
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